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ABSTRACT 

An up-to-date catalog of nearby galaxies considered as hosts of binary com- 
pact objects is provided with complete information about sky position, distance, 
extinction-corrected blue luminosity and error estimates. With our current un- 
derstanding of binary evolution, rates of formation and coalescence for binary 
compact objects scale with massive-star formation and hence the (extinction- 
corrected) blue luminosity of host galaxies. Coalescence events in binary compact 
objects are among the most promising gravitational-wave sources for ground- 
based gravitational- wave detectors such as LIGO. Our catalog and associated 
error estimates are important for the interpretation of analyses, carried out for 
LIGO, to constrain the rates of compact binary coalescence, given an astrophys- 
ical population model for the sources considered. We discuss how the notion of 
effective distance, created to account for the antenna pattern of a gravitational- 
wave detector, must be used in conjunction with our catalog. We note that the 
catalog provided can be used on other astronomical analysis of populations that 
scale with galaxy blue luminosity. 
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1. INTRODUCTION 

Compact binary coalescence (CBC) events, such as neutron star or black hole mergers, 
are one of the primary gravitational-wave sources for ground-based interferometers such as 
LIGOQ. LIGO's third (S3, Oct 31 2003 - Jan 9 2004) and fourth (S 4, Feb 22 2005 - M ar 23 



2005) science runs have reached significant extragalactic distances ( lAbbott et al.l 120071 ) into 
the nearby Universe. Especially for massive compact binaries whose components are black 
holes, the range extended beyond the Virgo Cluster. To interpret the searches for signals from 
compact binary coalescence in the LIGO data sets, it is necessary to use information about 
putative binary compact object populations in the known nearby galaxies, as well as how the 
population scales at larger distances. The nearby galaxy catalog discussed here represents 
the distribution of such extragalactic populations, and the procedures described are used for 
LIGO data analysis, such as assigning an astrophysically meaningful upper limits given non- 
detection. An accurate upper-limit that correctly incorporates our best information about 
galaxy distributions requires a model of the nearby overdense region because the current 
LIGO network's range probes this overdensity. 

Binary compact objects are usually produced from the evolution of massive stellar bi- 
naries. Since short-lived, massive stars emit more blue light than all other stars in a galaxy 
combined, blue light is a well-known tracer of star formation in general and the birthrate 
of these massive stars in particular. Given the short lifetimes of the known Milky Way 
double compact object populat ion and the slo w rate of change in star formation expected 



in nearby and distant galaxies, iPhinneyl (Il99ll ) has argued that a galaxy's blue luminosity 



should linearly scale with its compact binary coalescence rate. 

The sensitivity of LIGO to compact binary coalescence signals depends on the distance 
and sky position of the coalescence event and therefore, the distribution of known nearby 
galaxies in blue luminosity and in space is the minimum information needed to properly 
interpret searches of the LIGO data sets. 

It is possible that compact binary populations that are not related to regions of star 
formation may exist in the Universe. A mass, metallicity and morphology dependent star 
formation history may be needed to account for these population^. Nevertheless, the work 



1 www. ligo.org 



Lipunov et al.l (|l995f ) adopts a mass normalization to derive their event rate of 1/year within 50 Mpc 



using an older version of TuUy's catalog, whereas we use a blue-light normalization and the up-to-date TuUy 
catalog (ill). Their study also differs from ours because we consider issues like antenna pattern of the detector 
and completeness corrections, which they ignore. 
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described here is limited to the blue-hght luminosity as a tracer of the compact binary 
population. 

The contribution o f elliptical galaxies to the mer ger rates is potentially significant be- 
yond the Virgo cluster (Ide Freitas Pacheco et al.ll2006l ). whereas their blue luminosity is not 
representative of their putative compact binary populations. H owever, at large distan ces, the 
fractional blue luminosity produced in ellipticals is about 10% ( jPriver &: Allen! 120071 ). and at 
short distances the contribution is neglig i ble be cause there are fewer ellipticals in the nearby 
local universe. Ide Freitas Pacheco et al.l ( l2006l ) conclude that the event rate for an elliptical 
galaxy with the same blue luminosity as a spiral galaxy is a factor of five times larger on 
average. We conclude that LIGO rate upper limits derived from the catalog presented here 
would change by less than a factor of 1.5 due to a correction for elliptical galaxies. 

Our blue light census will also implicitly not account for any potential contribu tion from 
globul ar clusters to the compact binary coalescence rate of the nearby universe. iPhinney 
( ll99ll ) has argued that the contribution of globular clusters to double neutron star mergers 
in the Galaxy would not exceed 10 % of the coalescence rate due to the Galactic field. 
On the other hand it has been argued that the contri bution of globular clusters to binary 
black hole coalesceri ce may be very significant (see, e.g.. |Portegies Zwart and McMillarul2000l : 
O'Leary et al.ll2007l ). However, these cluster contributions are expected to become significant 
at distances beyond the Virgo cluster, where a more significant fraction of ellipticals with 
large globular-cluster systems will eventually enter LIGO's detection volume. 

We have used mostly publicly available astronomical catalogs of galaxies to compile 
a catalog used in the S3/S4/S5 (fifth science ru ifl, Nov 4 2005 - present) LIGO data set 
analyses. We discuss the methodology used to compile this galaxy catalog and briefly describe 
how this information feeds into LIGO rate estimates. In §21 we describe all the elements 
involved in compiling the galaxy catalog and assessing the relevant errors and uncertainties. 
In ^ we derive a correction factor to account for incompleteness in the catalog guided also 
by the blue-light volume density estimated from the Sloan Digital Sky Survey and earlier 
surveys. In §11 we discuss how the corrected catalog and resulting blue light distribution as 
a function of distance is used to bound the rate of compact binary coalescence using data 
from the recent LIGO science runs. If the maximum distance to which a search could detect 
a compact binary coalescence is known, then the expected number of detectable events can 
be derived. Some concluding remarks are made in §5] 
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COMPILATION OF GALAXY CATALOG 



We have compiled a cataloj^, the compact binary coalescence galaxy catalog or CBCG- 
catalog, of nearby galaxies which could host compact binary systems. For each galaxy out 
to 100 Mpc, the catalog provides the equatorial coordinates, distance to the galaxy, and the 
blue luminosity corrected for absorption. Estimates of the systematic errors on distance and 
luminosity are also provided. 

The CBCG-catalog is compiled from information provided in the following four astro- 
nomical catalogs: (i) t he Hubble Space Teles cope (HST) key project catalog used to measure 
the Hub ble constant (jFreedman et al.ll200ll ). (ii) Mateo's dwarf galaxi es of the local grou p 
catalog (lMatedll998l ). (iii) the HyperLeda (LEDA) database of gal axies (IPaturel et al.ll2003l ). 
and (iv) an updated version of the Tully Nearby Galaxy Catalog (iTuUyl l2006l ) . 



When combining these catalogs, distances and luminosities reported in the HST, Mateo 
and Tully catalogs were generally adopted over those in the LEDA catalog. This is because 
these catalogs use accurate distance determination methods compared to LEDA. Neverthe- 
less, LEDA served as the baseline for comparisons in the range 10-100 Mpc since it is the 
most complete. 



2.1. DISTANCES 

One of the primary objectives of the HST key project was to discover Cepheid variables 
(stars which have periodic variations in brightness) in several nearby spiral galaxies and mea- 
sure their distances accurately using the period-luminosity relation for Cepheids. Cepheid 
distance determination to nearby galaxies is one of the most important and accurate primary 
distance indicators. The distance information from the HST key project is considered to be 
the most accurate in the CBCG-catalog; there are 30 galaxies in our catalog for which we 
adopt distances from the HST key project. 



Mateo's review (lMatedll998l ) of properties of the dwarf galaxies in the Local Group pro- 
vides distance and luminosity information for each galaxy considered. Since the parameters 
in this catalog were derived from focused studies on each individual galaxy, we consider it 
the most accurate next to the HST measurements for nearby galaxies. Moreover it has rea- 
sonably comprehensive information on the Local Group's dwarf galaxies; there are 18 sources 



^http://ww w. lsc-group.phys.uwm.edu/cgi-bin/cvs/viewcvs.cgi/lalapps/src/ 
inspiral / inspsrcs 1 OOMpc. errors?cvsroot=lscsoft 
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in the CBCG-catalog which adopt distances (and luminosities) from Mateo's compilation. 

It becomes increasingly difficult to use primary distance estimators like Cepheid stars 
in more distant galaxies. Therefore secondary distance methods are used to measure larger 
distances. Tully's catalog has up to three types of distances for each source: (i) Quality 
distance (-Dq) is based on either Cepheid measurements, surface brightness fluctuations, or 
the tip of the red giant branch. There are 409 galaxies with such a distance in the CBCG- 
catalog. (ii) HI luminosity-line-width distances (-Dhi) are obtained from the TuUy-Fisher 
relation, where the maximum rotational velocity of a galaxy (measured by the Doppler 
broadening of the 21-cm radio emission line of neutral hydrogen) is correlated with the 
luminosity (in B, R, I and H bands) to find the distances. There are 553 galaxies in the 
catalog with such a distance, (iii) Model distance (-Dm) is derived from an evolved dynamical 
mass model that translates galaxy r adial velocities into distances. This model is an update of 



the least action model described by lShaya et al.l (119951 ) and takes into account the deviations 



from a perfect Hubble flow due to a spherically symmetric distribution of mass centered on 
the Virgo Cluster. All galaxies have a calculated model distance. Whenever available, 
distances are the most preferred due to their smaller uncertainties, then the Dm followed by 
Dm. 

The remaining galaxies come from LEDA which does not provide distances explicitly, 
but instead provides measured radial velocities corrected for in-fall of the Local Group to- 
wards the Virgo cluster (fvir)- We obtain the LEDA dista nce {D-[) using Bubb le's law with 



the Hubble constant Hq = 73kms~^Mpc~^ reported by ISpergel et al.l (120061 ). Although 
corrections to the recessional velocity were made, this method of calculating distances is 
still highly uncertain. Hence, we use Bubble's law to evaluate the distances only to the 
galaxies for which v^i,- > 500 km/s (7Mpc) and peculiar velocities are expected to be more 
of a perturbation. 

The error in a distance depends strongly on the method used to measure that distance. 
The HST sources, though a sma ll contribution to the galaxy catalog, have the smallest errors 
(< 10%) ( Freedman et ahlboOl ). The three different distance methods in Tully's catalog have 



different errors. Dq also has a low error (10%) followed by the Dm (20%). To obtain an 
estimate for the errors of Dm, we compare them with Dq for the set of galaxies that have 
both types of distance estimates. The best fit Gaussian (see Fig. [1]) to the logarithm of 
fractional errors has a one sigma width of 0.24 which when subtracted in quadrature with 
Dq error gives, 0.22 distance error associated with Dm- 

Because errors in v^i,- are not given in LEDA, we follow a similar procedure to find LEDA 
distance errors, Dl- We compare the calculated Dl with Dq for galaxies in both catalogs 
to obtain uncertainty estimates in Dl. The plot in Fig. [2] shows the best fit Gaussian to the 
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logarithm of fractional errors with a one sigma width of 0.27 which, subtracted in quadrature 
with Dq distance errors, gives a total distance error 0.250 



^For searches of the S3 and S4 LIGO data (jAbbott et al. 
uncertainty of 40% was used for LEDA distances. 



20071 ) ■ with smaller ranges a more conservative 
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Fig. 1. — In order to obtain reasonable estimates for Tully's model distances we compare 
galaxies that have values for both. We only consider galaxies beyond 10 Mpc since model 
distances and LEDA distances are not reliable below this value. All galaxies below 10 Mpc 
have better distance estimates. The Tully quality distance has roughly a 0.1 logarithmic 
error. The best fit Gaussian for \n\DM / Dq\ implies a fractional error a of 0.24 in log. 
Subtracting these uncertainties in quadrature gives an error of 0.22 for Tully model distances. 
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Fig. 2. — Fractional error analysis as in Fig. [T] for LEDA distances. By comparing the 
fractional error between LEDA distances and Tully we obtain a ~ 0.25 log distance error for 
LEDA. 
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2.2. BLUE LUMINOSITIES 



The distribution of binary compact objects in the nearby universe is expected to follow 
the star formation in the universe and a measure of star forniation is the blue luminosity of 
galaxies corrected for dust extinction and reddening (jPhinneylll99ll ). Hence, for each galaxy, 
we calculate the blue luminosity Lb from the absolute blue magnitude of the galaxy Mb 
(corrected for internal and Galactic extinctions). For convenience, blue luminosity is provided 
in units of Liq = IO^^Lb,©, where Lb,q = 2.16 x 10^^ ergs/s is the blue so lar luminosity derived 
from the blue solar magnitude Mb,© = 5.48 (IBinney fc Tremainell2000l ). We do not consider 
galaxies with luminosities less than 10"'^Lio because they do not contribute significantly to 
the total luminosity - see ^ 

The Mateo, Tully and LEDA catalogs provide information on apparent B-magnitudes 
corrected for extinction. The galaxies in the HST key project catalog have only distance 
information, so for those we extract the corresponding apparent magnitude values (ttib, 
corrected for internal and Galactic extinction) in the B-band from the Tully catalog to find 
Mb. Table [T] summarizes relevant properties of each of these catalogs and the fraction of the 
total luminosity within 100 Mpc that each contributes. 
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Table 1: Summary information about the four astronomical catalogs used to develop the 
CBCG-catalog. We report the number of galaxies for which the catalog was the primary 
reference and fraction of the total CBCG-catalog blue luminosity accounted for by those 
galaxies. 



Catalog 


# of galaxies 


Lio 


Fractional luminosity 


Reference 












(i) HST 


30 


57.3 


0.1% 


fFreedman et al. 2001) 


(ii) Mateo 


18 


0.4 


<0.001% 


fMateo 1998) 


(iii) Tully 


1968 


2390 


5.3% 


fTuUv 2006) 


(iv) LEDA 


36741 


42969.4 


94.6% 


fPaturel et al. 2003) 


Total 


38757 


45417.1 


100.0% 
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The LEDA database quotes uncertainties in apparent magnitude. Figure [3] shows the 
distribution of LEDA assigned apparent magnitude variances for the galaxies in the CBCG- 
catalog. The RMS error is Am- R = 0-42. Galaxies from Tully's catalog have a smaller 
observational error Ams = 0.30 (lTullyll2006l ). 



12000 
10000 




Fig. 3. — LEDA provides uncertainties in apparent magnitudes. The histogram above 
the 777.6 variance distribution for each LEDA galaxy. The RMS error is 0.42. 
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COMPLETENESS 



Observations of faint galaxies are difficult even in the nearby universe and lead to system- 
atic incompleteness in galaxy catalogs. Studies of galaxy luminosity functions can provide 
insight into how many galaxies are missing from a catalog (and hence the corresponding 
blue luminosity). Using the CBCG-catalog, we can generate a luminosity function N{L,D) 
which is the number of galaxies with luminosities within a luminosity bin from L to L + AL 
normalized to the spherical volume within radius D. Specifically, we write 



N{L,D)AL 



where 



if (L < Lj < L + AL) and {Dj < D) 
otherwise 



and the sum over j runs through all the galaxies in the catalog. The quantities Lj and D-^ 
are the luminosity and distance of each galaxy. Similarly we can compute the luminosity 
function in terms of blue absolute magnitudes as a function of distance N{Mb, D). The 
dashed and dot-dashed lines in Fig. H] show several realizations of N{Mb,D) for different 
distances D plotted as a function of Mb. 

To estimate the degree of incom pleteness in the CBCG-catalog, we use an analytical 
Schechter galaxy luminosity function (jSchechterl 119761 ) 



(2) 



where 0(L)(iL is the number density (number of galaxies per unit volume) within the lumi- 
nosity interval L and L + dL, L* is the luminosity at which the number of galaxies begins 
to fall off exponentially, a is a parameter which determines the slope at the faint end of 
the luminosity function, and (f)* is a normalization constant. In terms of (blue) absolute 
magnitudes, Mb, the Schechter function becomes 



4>{MB)dMB = O.92 0*exp [-iQ-^'-^i^'^^-M*)^ ^^q- 



-0.4{MB-Af*) 



a+1 



dM 



B • 



(3) 



To estimate the total lumi nosity function , we u se results from the Sloan Digital Sky 
Survey (SDSS) as reported by iBlanton et al.l (120031 ). Although the SDSS sky coverage 
is inadequate in RA and DEC, it provides excellent coverage throughout our desired dis- 
tance and beyond. We therefore use the green luminosity function Schechter fit given 
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in Table 2. of iBlanton et al.l ( 20031) and conve rt it into blue band using the expression 



given in Table 2. of iBlanton &: Roweij (120071 ). Adopting a Hubble constant value of 



73kms ^ Mpc ^ ( Spergel et al. 2006 ) and correcting for reddening]^ the Schechter param- 
eters are {M^,(p*,a) = (-20.3,0.0081,-0.9). The solid line in Fig. H shows the Schechter 
function 0(Mb) derived from these values. Since this function is obtained from deep sur- 
veys, it does not account for the local over-density of blue light coming primarily from the 
Virgo cluster. For distances up about to 30 Mpc, the CBCG-catalog's luminosity function 
N{Mb,D) exceeds 0(Mb). 

We can now derive a completeness correction that arises at the faint end beyond about 30 
Mpc, where the Schechter function exceeds the catalog N{Mb, D). We integrate the CBCG- 
galaxy-catalog luminosity function N{L, D) over L and subtract it from the Schechter fit 
as a function of distance. Hence, the total corrected cumulative luminosity Ltotai within a 
volume of radius D is given by 

LtoUD) = Lcbcg(^) + L,,rr{D) (4) 

where 

Lcbcg{D) = / dD'Y^LAD'-D^) (5) 

Jo : 



^max 



Air f^" 

L,,rr{D) = —D' I LdLQ[(l){L)-N{L,D)] [(l){L)-N{L,D)] . (6) 

-^min 



3 



Here, the index j runs through all galaxies in the catalog, 6 is the Dirac delta function, G is 
the step function and 4>{L) is the adopted Schechter function (distance independent) assumed 
to represent the complete luminosity distribution. We note that Lmax = 52.481 Lio (Mb = 
—23.83) is the maximum luminosity in the CBCG-catalog and we choose L^i^ = 10~^Lio 
(Mb = —12.98) because luminosities below this value do not contribute significantly to 
the net luminosity. The quantity Lq-qcg in Eqs. and (jSD is the uncorrected cumulative 
luminosity from the CBCG-catalog; the quantity Lcott is the completeness correction. Note 
that the completeness correction term is always zero or positive regardless of the choice of 
Schechter function. 

In Fig. O we show the cumulative blue luminosity as a function of distance as obtained 
directly from the CBCG-catalog (solid line) as well as with the completeness correction 
applied (dashed line). It is evident that the correction becomes significant at distances in 
excess of about 40Mpc. 



^We correct the value of Afg to be consistent with the reddening correction described in H3.ll 
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Fig. 4. — The luminosity function of CBCG catalog at various distances (d ashed and dot- 



dashe d lines) and a Schechter function fit (solid line) given in Eq. ([3]) based on lBlanton et al. 



(120031 ) ■ We compensate for the incompleteness of the CBCG-catalog by applying an upward 
correction to the luminosity bins that are below the Schechter function fit (solid line), ac- 
cording to Eqs. (jl]) and ([6]). Error bars are found by sliding the magnitudes of each galaxy 
according to the mean errors and recomputing the luminosity function. 
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3.1. Comparison with other results 



To compare our method of correcting for completeness with other metho ds, we consider 



the di rect computation of a reddening corrected luminosity density based on iBlanton et al. 
(I2OO3I ) which could be used at large distances. We adopt a blue luminosity density of 
(1.98 ± 0.16) X lO-^Lio/Mpc^ calculated as follows: 



The blue luminosity density, in terms of blue absolute magnit udes per cubic Mpc, is 
-14.98 locally (redshift 2 = ) and -15.17 for z = 0.1 [Table lo lsianton et al.l J2003h ]. 
This is for a standard cosmology with Qu = 0.3 and I^a = 0.7. We use z = 0.1 so that 
the results will be valid for advanced detectors. 

We convert the z = 0.1 blue magnitude density (-15.17) to luminosity units 1.33 x 
10~^Lio/Mpc^ and assign systematic errors (~ 10%) associated with the photometry 
to obtain a luminosity density of (1.33 ± 0.13) x 10~^Lio/Mpc^. 

We also correct for pr ocessing of blue li ght an d re-emission in the infrare d (IR) following 
PhinnevI (Il99ll ) and iKalogera et al.l (1200 ll ) . We use the analysis of ISaunders et al. 
(119901 ). upward correct by 30% th eir far IR (AOp^m — l OO/xm) luminosity density to 
account for emission down to 12/im (IKalogera et al.ll200ll ). and convert to Liq to obtain 
an IR luminosity density of Ljr = (0.65 ±0.1) x 10~^Lio/Mpc^. 

Adding both luminosity densities above and accounting for the errors, we obtain a blue 
light luminosity density corrected for extinction equal to (1.98 ±0.16) x lO^^Lio/Mpc^ 



We use this blue luminosity density and its uncertainty and plot the implied cumulative blue 
luminosity as a function of distance (cubic dependence) in Fig. [5] (gray-shaded region). This 
uniform density distribution agrees well with the completeness corrected luminosity given 
above. 

We can compare our resu lts for the cumulative b lue luminosity as a function of distance 
to similar results obtained by iNutzman et al.l (120041 ). especially their Figure 1. The results 
for the uncorrected catalog agree qualitative ly. However, the catal og described here is more 
up-to-date compared to the one compiled by iNutzman et al.l (120041 ) by virtue of the updates 
to LEDA and by the inclusion of the current Tully catalog. The incompleteness correction 
derived here is also more physically and empirically motivated than the one constructed in 
that earlier paper. We note that the cumulative luminosity shown as the dashed line in 
their Figure 1 is too low by a factor of 47r/3 due to a numerical error. Additionally, their 



luminosity d ensity is ~ 25% 



presented by iBlanton et al.l (120031 ) 



ower than ours resulting from our use of the more recent results 
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Distance (Mpc) 

Fig. 5. — Cumulative luminosity as a function of distance from CBCG-catalog uncorrected 
for incompleteness (solid line), corrected for incompleteness (dashed line) and the cubic 
extrapolation from the assumed constant blue luminosity density corrected for extinction 
(gray-shaded region). 
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COMPACT BINARY COALESCENCE RATE ESTIMATES 



For neutron star binaries, the observed binary pulsar saniple can be used to predict 
the coalescence rate TZmw in the Milky Way (IKim et al.l 12004 |2006| ) . The coalescence rate 
within a sphere of radius D is then simply given by 



R = n 



MW 



^total 



(7) 



where Ltotai(-D) is the total blue l uminosity within a d istance D and Lmw is the blue lu- 
minosity of the Milky way, I.TLiq (IKalogera et al.ll200ll ). If the rate i? of a binary neutron 
star coalescence could be measured directly, it would provide an independent estimate of 
the rate of coalescence per unit of blue luminosity. Together these two measurements would 
deepen our understanding of stellar and binary evolution. Furthermore, the current un- 
derstanding of binary evolution and compact object formation lead s us to anticipate the 



form a tion of black hole binar ies that will merge within a Hubble time (je.g.. Belczynski et al. 



2OO2I : iBelczynski et al.l 120071 ). Experiments like LIGO will provide a direct measure of the 
compact binary coalescence rate and will impose constraints on the theoretical models of 
stellar evolution and compact binary formation. 



4.1. Rate estimates and systematic errors in gravitational- wave searches 

In its simplest form, the rate estimate derived from a gravitational- wave experiment will 
take the form 

constant , , 

= 7^77^ — (8) 

TCl ^ ' 

where the constant depends on the precise outcome of the search and the statistical method 
used in arriving at the rate estimate, Cl is the cumulative blue luminosity observable within 
the search's sensitivity volume measured in Lio, and T is the time analyzed in years. In 
general the sensitivity volume is a complicated function which depends on the instrument 
and the gravitational waveforms searched for. Here, we focus on the influence of the host 
galaxy properties and the distribution of blue light with distance. 

The gravitational-wave signal from a compact binary inspiral depends on a large number 
of parameters. It is convenient to split these parameters into two types for our discussion. 
Of particular interest here are the parameters which determine the location and orientation 
of the binary. We denote these collectively as A := {D,a,5, L,ip,t}, that is the distance 
to the binary, its Right Ascension and declination, inclination angle relative to the line of 
sight, polarization angle of the waves, and the time when the binary is observed, respectively. 
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Other parameters, including the masses and the spins, are denoted jl. Recognizing that the 
spatial luminosity distribution can be written as 



L{a, 6,D) = J2 S{aj - a) 6{6j - 6) 6{Dj - D) 



(9) 



we write the cumulative luminosity as 



Cr 



L{a, 6, D) p(detection|/I, A) p{l) pijp) pit) djldX 



(10) 



Assuming that binary coalescences are uniformly distributed in time, and their orientation 
is random, we take the corresponding prior probabilities: p{l) = sin(i)/2, p{t) = 1/day, and 
piip) = l/27r. 

Systematic errors associated with the derived rate esimates are naturally associated 
with the errors in cumulative luminosity Cl- The two most relevant errors in the galaxy 
catalog are in apparent magnitude rriB and distance D. Sky positions are known so precisely 
that small errors in RA and DEC do not change the detection probability of a particular 
binary in any significant way ; for this reason, such errors are not included in the LIGO 
analyses (lAbbott et al.l 120071). The errors induced on the spatial luminosity function in 
Eq. ([9]) take the form JPairhurst et aliboOTh 



AD 



[L+AL]{a,6,D) = io-°-^^"^s, 1 1 + __i J S(aj~a) 6{5j-6) 6{Dj+ADj-D) . (11) 



4.2. A simplified model for estimating expected event rates 

The sensitivity of a search for gravitational waves from compact binary coalescence is 
determined primarily by the amplitude of the waves at the detector. For a non-spinning 
binary (i.e., the spins of each compact object are much smaller than their general- relativistic 
maximum value of mf) with given jj, the a mplitude is inversely proportional to the effective 
distance D^^ defined as ( lAUen et al.l 120051 ) 



D 



cff 



V/F2(1 + C0S2 6)2/4 + F2C0S2 



(12) 



where D is the physical distance to the binary, F+ and Fx are the response amplitudes 
of each polarization at the detector which depend upon the location of the binary system 
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flAnderson et al.ll200l[ ): 

= ~2 ~'~ "^"^^^ cos 20 cos 2'?/' — cos 6 cos 20 sin 2?/' (13) 
Fx = ^(l + cos^6')cos20sin2^/'-cos6'sin20cos2^/^ . (14) 

Here 9 and are the spherical co-ordinates of the source defined with respect to the detector 
and, as before, l and ip are the inchnation and polarization angles. Since 6 and are detector 
dependent, the effective distance is different for geographically separated detectors that are 
not perfectly aligned and, for a fixed source location, changes as the Earth rotates through 
a sidereal day. Additionally, the effective distance is always at least as large as the physical 
distance. 

For simplicity in understanding the sensitivity of gravitational-wave searches, consider 
the case in which /lis fixed, i.e. j9(/I) = 5{Jl—jl). For example, these might be the parameters 
appropriate to a neutron star binary. The sensitivity of a detector is given by its horizon 
distance, which is defined as the maximum effective distance that a neutron star binary 
system can be detected at signal-to-noise ratio of 8. Consider a search which can perfectly 
detect these binaries if they have an effective distance D^s < -Dhorizon at a particular detector. 
Then 

— * — * 

p(detection|/i. A) = 6(Dcfr(A) < F'horizon) (15) 

and we can write 

CL(/^horizon) = j Lia,5,D) e(Deff(A) < Dhorizon) P(Op(^)p(^)^^A. (16) 

Thus, the cumulative blue luminosity accessible to such a detector is the blue luminosity 
within an effective distance sphere of radius -Dhonzon, averaged over the time of day and 
possible orientations of the binary. The lower curve in Fig. |6] shows C/, (/^horizon)- Figure 
E] also illustrates the significant difference between the cumulative luminosity Cl{Dcs) and 
total luminosity Ltotai(-D) at a given distance. If galaxies are distributed uniformly in space 
the ratio between these is ~ 11.2 ; this is the factor by which the detection rate would be 
reduced and arises purely from the LIGO detector response, averaged over all possible source 
orientations with respect to the detector. 

When estimat ing the rate based on gravitational- wave observations, one can marginalize 



over uncertainties (IFairhurst et al.ll2007l ) in the galaxies' distances and apparent magnitudes. 
Specifically, by making use of the modified spatial distribution function Eq. (fTT!) and the 
distributions for ADj and Amsj reported here, we can obtain a probability distribution for 
the cumulative luminosity p(Ci| Am^j, ADj) from Eq. fllOp . For each value of the cumulative 
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luminosity, a probability distribution p{R\Cl) for the event rate can be calculated. Finally, 
the rate is marginalized over errors in the galaxy catalog by computing 



This distribution is then used to obtain a rate interval or upper limit on the occurrence of 
binary coalescences in the unverse. 

While this approach provides a reasonable estimate of the observable blue light luminos- 
ity in a single detector, it does not provide the whole story. For example, the 16° difference in 
latitude between the LIGO Observatories in Hanford, Washington and Livingston, Louisiana, 
implies the Ci(Dhorizon) depends on the site used. Figure [7| shows two-dimensional contours 
of this function. 




(17) 



-22- 




lo" 10^ 10^ 10^ 

Distance (Mpc) 



Fig. 6. — Cumulative luminosity as a function of physical distance (top line) and horizon 
distance (bottom line). The horizon distance -Dhorizon is defined as the physical distance to 
an optimally oriented and located binary system that would be detected with a signal-to- 
noise ratio of 8. (Instrumental sensitivity range is sometimes quoted in terms of the radius 
of a sphere with the same volume as the non-uniform region probed by the instrument, 
this sensitivity range Dg is related to the horizon distance by Dg ~ -Dhorizon/ "\/5- The gray 
shaded lines are cubic extrapolations (^ derived for both cases. Given a LIGO horizon 
distance one can immediately get the cumulative blue luminosity from the bottom curve. To 
obtain an approximate rate upper limit one could calculate 7^9o% [j'^^^L'^^] = 2.3/(Ci x T) 
where Cl is taken from this plot at a given range in horizon distance. Inset: Ratio of the 
cumulative luminosity for the physical and horizon distance from the completeness corrected 
CBCG-catalog illustrates the non-uniform distribution at smaller ranges (< 20 Mpc) and 
asymptotes to the expected uniform distribution ratio (dashed line) for larger distances. 
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Fig. 7. — Luminosity contours per effective distance bin in the two LIGO sites. The effective 
distance to a source in one galaxy is different between tlie two detectors, cliang 
function of the sidereal day and also on the orientation of the particular source. Since the 
effective distance is always larger than the real distance the luminosity available within a 
given effective distance bin is considerably smaller than the luminosity within the physical 
distance bin. The upper horizontal numbers refer to the luminosity per bin in effective 
distance. The parenthetical lower numbers refer to the luminosity per physical distance 
bin. It is also possible to have a systematically different luminosity between detectors as 
is indicated in the right panel zoom of the first 5 Mpc. The available luminosity within 5 
Mpc (mostly from Andromeda) is slightly better located for LLO and therefore stretches 
the contours to higher effective distances for LHO. LIGO rate upper limits for searches with 
limited range thus depend on the non-uniformity of the Local Group. 
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Based on the galaxy catalog presented in this article, the cumulative blue luminosity 
Cl, measured in Lio, accessible to a search with a given horizon distance sensitivity can be 
derived from Fig. O and is tabulated in Table [2J We can combine the calculated cumulative 
blue luminosity with estimates of TZ, the rate of binary mergers per Lio, to estimate the 
number of compact binary merger events detectable in a given LIGO search with an 
observation time T: 



N =10 



-3 



7^ 



Ci 



10 



T 

yr 



If the horizon distance of a search is larger than 50 Mpc, we can use the following approxi- 
mation, from a cubic law: 



N ^7Ax IQ- 



n 



^10 



Myr" 



-^horizon 

100 Mpc 



)■(!) 
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Estimated rates of binary neutron star (BNS) mergers in our Galaxy are based on 
the observed sample of binary puls a rs. T he rates depend on the Galactic distribution of 
compact objects. In iKalogera et al.l (120041 ). the most recent reference estimating rates, the 
most likely Galactic rate for their reference model 6 is 83Myr~^, with a 95% confidence 
i nterv al 17 — 292Myr~^. The most likely rates for all the models used in IKalogera et al. 
are in the range 4 - 220 Myr"^ for the Milky WayQ 



For the 4km LIGO detectors currently operating, /^horizon ~ 30 Mpc for BNS. Thus, the 
predicted number of BNS events is in the range iVg ~ 2 — 30 x 10~^ yr~^ with the most likely 
number being iVfi ^ 1/(100 yr ) [we u se the subscript g to indicate these rates use reference 
model 6 from IKalogera et al.l (|2004l )]. A search that reaches twice the distance (such as 
enhanced LIGO), yields a most likely rate 1/(10 yr). And a search that would be 

15 times more sensitive to coalescences of binary systems than the current LIGO detectors 
(such as Advanced LIGO) would yield a most likely rate of Nq ^ 40.0 yr^^. 



'''The rates quoted here are in units of rate per Milky Way per Myr; to get the rate per Lio, we divide 
by 1.7 which is the estimated blue lumino sity of the Milky way in Lio units, assuming the blue absolute 
magnitude of the Milky Way to be -20.11 |Kalogera et al.lboOll ). 
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Table 2: Table showing the cumulative blue luminosity Ci(-Dhorizon) accessible to a search 
with horizon distance -Dhorizon given in the first column. For -Dhorizon > 100 Mpc, the cu- 
mulative blue luminosity accessible to a search is given approximately by Cl (Dhorizon) ~ 

7.4 X 10^ (L»horizon/l00Mpc)3. 



-Dhorizon (Mpc) 


Cl( Dhorizon) (-^lo) 


10 


23 


20 


85 


30 


240 


50 


953 


100 


7200 


200 


59200 


300 


200000 


500 


926000 
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5. CONCLUSION 

Whether one wishes to compute expected detection rates for LIGO searches, or to 
interpret LIGO searches as rate upper hmits (or eventually detection rates), we require at the 
simplest level accurate accounting of the total observable blue luminosity Cl. As mentioned in 
the previous sections, a galaxy catalog complete with sky positions and distances is important 
for first generation LIGO detectors because the blue luminosity is not uniformly distributed 
in the sky within the search range. An upper limit which takes in to account the most up- 
to-date information on galaxy distribution can be obtained by accurately modeling the local 
overdense region. For searches with ranges well beyond current sensitivity the universe is 
uniform and rate estimates depend primarily on accurate blue luminosity densities corrected 
for reddening. We have introduced a method to bridge the gap between the well known 
nearby galaxy distribution and the expected lon g range distribution through a completeness 



correction based on SDSS luminosity functions (IBlanton et al.ll2003l ). 



This paper provides the most up to date accounting of nearby galaxies within lOOMpc 
as well as errors in the apparent magnitude (corrected for reddening) and distance and 
demonstrates how the errors propagate into rate calculations. Astrophysical errors are a 
significant co ntribution to the eve ntual systematic error associated with coalescence rate 



upper limits (IFairhurst et al.l 120071 ) and must be included. This paper provides a survey of 
the asymptotic and local uncertainty. Motivated by the use of effective distance to account for 
the antenna pattern of the LIGO detectors, we demonstrate the need to compute the average 
blue light luminosity within a given effective distance sphere. For ranges within 50Mpc there 
is a nontrivial relationship between cumulative blue luminosity within an effective distance 
sphere and within a physical distance sphere. Beyond 50Mpc the relationship is well behaved 
leading to the simple scaling for the number of detected events N given in Eg. (1191) . We would 
like to point that the catalog provided can also be used on other astronomical analysis of 
populations that scale with galaxy blue luminosity, such as the local Type II supernova rate 
or the rate of nearby SGR bursts that show up as short GRBs. 

We provide sufficient description of our methods for others to apply new rate models 
to future LIGO data. Although this catalog will serve as a reference for current and future 
LIGO data analysis, we look forward to future work that may transcend the simple blue 
light rate normalization that we have discussed. One way to go beyond blue light rate 
normalization, (necessary to ascertain the degree to which old stars contribute to present 
day mergers) is with multiband photometry of nearby galaxies which can reconstruct their 
mass, morphology and metallicity dependent star formation history. With this information 
in hand LIGO detections could be applied more stringently to assess the relative contribution 
that progenitors of different ages provide to the present day merger rate. 
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